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09/407,617, filed on September 28, 1999; 09/163,933, filed on September 29, 1999; and 
09/461,787, filed on December 15, 1999, the disclosures of which are hereby 
incorporated by reference as if fully set forth herein. 



Field of the Invention 

The invention described below generally relates to controlling and diagnosing the 
health of a machine, and more particularly, to systems and methods for controlling and 
diagnosing motorized systems according to vibration, pressure, temperature, speed, 
1 5 and/or current analysis. 



Background of the Invention 

Many industrial processes and machines are controlled and/or powered by electric 
motors. Motorized systems include pumps providing fluid transport for chemical and 

20 other processes, fans, conveyor systems, compressors, gear boxes, motion control 

devices, screw pumps, and mixers, as well as hydraulic and pneumatic machines driven 
by motors. Such motors are combined with other system components, such as valves, 
pumps, conveyor rollers, fans, compressors, gearboxes, and the like, as well as with 
appropriate motor drives, to form industrial machines and actuators. For example, an 

25 electric motor may be combined with a motor drive providing electrical power to the 
motor, as well as with a pump, whereby the motor rotates the pump shaft to create a 
controllable pumping system. 

Controls within such motorized systems provide for automatic system operation 
in accordance with a setpoint value, and may additionally allow for manual operation. 

30 Thus, for instance, a motorized pump system may be operated so as to achieve a user 

specified outlet fluid flow rate, pressure, or other system setpoint. In another example, a 
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motorized conveyor system may include one or more motorized roller systems, wherein 
the individual roller systems are controlled according to a conveyor speed setpoint Such 
motorized system controls may include a controller receiving a setpoint from a user or 
from another system, which inputs one or more system performance values and provides 
5 appropriate control signals to cause the motorized system to operate in a controlled 
fashion according to a control scheme. For example, a motorized pump system may be 
controlled about a flow rate setpoint, wherein the associated controller reads the setpoint 
from a user interface, measures the system outlet flow rate via a flow sensor, and 
provides a control signal indicative of pump speed to a motor drive operatively connected 

1 0 to a motorized pump, whereby the control signal is adjusted so as to achieve the setpoint 
flow rate in closed-loop fashion. 

Although operation of such motorized systems and controllers may achieve 
system operation in accordance with the setpoint, other factors such as system component 
wear, component faults, or other adverse conditions, and the like, may affect the 

15 operation of the motorized system. Thus, for example, degradation in a pump impeller in 
a motorized pumping system may lead to premature catastrophic failure of the system if 
left unchecked. In this regard, operation of the pump strictly in accordance with a flow 
rate setpoint may accelerate the system component wear, degradation, and/or failure, 
whereas operation at other flow rates may allow the system to last longer. This may be of 

20 importance in critical systems where safety is an issue. For instance, the motorized 
pumping system may be located on board a military vessel at sea, wherein operation 
according to a flow setpoint may lead to catastrophic pump impeller failure before the 
vessel can be brought to port for repairs or maintenance, whereas system operation at a 
reduced flow rate may allow the pump to survive until the next scheduled servicing. 

25 Motor diagnostics apparatus has been employed in the past to provide an 

indication of wear, damage, and/or degradation in motorized system components, such as 
motors, prior to catastrophic failure thereof. Such diagnostic devices may be used to 
monitor the overall health of either the motorized system components being controlled, or 
the control system itself. In this regard, assessing system health can be used to minimize 



-2- 



01AB121 



unscheduled system downtime and to prevent equipment failure. This capability can 
avoid a potentially dangerous situation caused by the unexpected outage or catastrophic 
failure of machinery* However, many conventional diagnostic devices inconveniently 
require an operator to manually collect data from machinery using portable, hand-held 
5 data acquisition probes . 

Other known systems have sensors and data acquisition and network equipment 
permanently attached to critical machinery for remote diagnostics. Typically the 
diagnostics equipment is directed to detecting problems with the control system hardware 
itself or monitoring the integrity of the output, i.e., monitoring when the control system 

10 response is outside prescribed time or value limits. As noted above, control system 

health monitoring, health assessment and prognostics generally are performed in isolation 
from any associated control system. These systems typically conduct passive monitoring 
and assess system health using diagnostic algorithms and sensors dedicated to establish 
system health. This passive monitoring is frequently done using off-line, batch-mode 

1 5 data acquisition and analysis to establish the health of the system. 

In conventional motorized systems, therefore, controlled operation is provided 
about a setpoint, wherein such controlled setpoint operation may exacerbate system 
component degradation and/or accelerate catastrophic failure thereof. Prior diagnostics 
apparatus achieves some level of identification of such system component degradation 

20 prior to component failure. However, as noted previously, because virtually all 
diagnostics systems perform off-line diagnostic processing, it has been extremely 
difficult to implement diagnostics processing real-time in coordination with on-line 
control. Thus there is a need for improved control and diagnostics systems and 
techniques by which controlled operation of motorized systems can be achieved while 

25 mitigating the extent of component degradation and failure. 



Summary of the Invention 

The following presents a simplified summary of the invention in order to provide 
a basic understanding of one or more aspects of the invention. This summary is not an 
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extensive overview of the invention. It is intended to neither identify key or critical 
elements of the invention, nor to delineate the scope of the present invention. Rather, the 
sole purpose of this summary is to present some concepts of the invention in a simplified 
form as a prelude to the more detailed description that is presented hereinafter. The 
5 invention provides control systems and methodologies for controlling and diagnosing the 
health of a motorized system and/or components thereof. Diagnosis of the system or 
component health is accomplished using advanced analytical techniques such as neural 
networks, expert systems, data fusion, spectral analysis, and the like, wherein one or 
more faults or adverse conditions associated with the system may be detected and/or 
10 diagnosed. 

The diagnostics may be performed using one or more measured system 
parameters, such as vibration, pressure, temperature, speed, power, current, or the like. 
31 Frequency spectral analysis may also be employed in order to detect and/or diagnose 

component wear, degradation, failure, faults, and the like. For example, one or more 
15 signals from system sensors may be processed by a diagnostics system and analyzed in 
the frequency domain, whereby such adverse conditions may be identified and/or 
predicted. The diagnostics system may further provide one or more diagnostics signals 
indicative of the health of a motorized system, which may then be provided to an 
associated controller, whereby the operation of the motorized system may be modified. 
20 In this regard, the system operation may be modified to ameliorate detected component or 
system degradation. 

The invention may thereby provide for extending the useful service life of one or 
more system components. Thus, the invention finds utility in association with a variety 
of industrial processes and machines which are controlled and/or powered by electric and 
25 other types of motors. Such processes and machines include pumps providing fluid 
transport for chemical and other processes, fans, conveyor systems, compressors, gear 
boxes, motion control devices, screw pumps, and mixers, as well as hydraulic and 
pneumatic machines driven by motors. 

According to an aspect of the present invention, there is provided a diagnostics 
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and control system for controlling a motorized system and diagnosing the health thereof. 
The diagnostics and control system may comprise a controller operatively associated with 
the motorized system and adapted to operate the motorized system in a controlled 
fashion, and a diagnostics system operatively associated with the motorized system and 
5 adapted to diagnose the health of the motorized system according to a measured attribute 
associated with the motorized system. The measured attribute may comprise vibration, 
pressure, current, speed, and/or temperature, for example, wherein such attribute 
information is obtained from one or more sensors operatively connected to the motorized 
system. 

1 0 The diagnostics system may provide a diagnostics signal according to the health 

of the motorized system, and the controller may accordingly provide a control signal to 
the motorized system according to a setpoint and/or the diagnostics signal. For example, 
the controller may comprise one or more changeable parameters, such as gains or the 
like, wherein the parameter may be changeable in response to the diagnostics signal. In 

1 5 this manner, the operation of the motorized system may be adaptable to a variety of 

system health conditions, whereby the adverse effects of system health problems may be 
proactively addressed in order to mitigate the effects thereof, including for example, the 
extension of service life and the reduction in system downtime and/or failures. 
Moreover, the output of the controller may be provided to the diagnostics system so that 

20 the health assessment made by the diagnostics system can be based at least in part on the 
control signal and the response by the motorized system to this control action. 

In addition to component diagnostic information, the invention provides for 
detecting, diagnosing, and/or adaptive control according to other detected system 
operating conditions. For instance, the diagnostics system may be advantageously 

25 employed in order to detect and diagnose cavitation, blockage, and/or the like in a 
motorized pump or in the system in which the motorized pump is employed. In this 
regard, pressure and flow information related to the pumping system may be measured 
and cavitation conditions may be identified using a classifier system, such as a neural 
network. The diagnostics system may thus comprise such a classifier system for 
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detecting pump cavitation according to flow and pressure data. The invention may be 
employed in cavitation monitoring, as well as in control equipment associated with 
pumping systems, whereby pump wear and failure associated with cavitation conditions 
may be reduced or eliminated. 
5 Such pumping system cavitation conditions may alternatively or in combination 

be diagnosed using other measured signals from the system. For instance, it has been 
found that fault data relating to the operating condition of a motorized pump may be 
ascertained from variations in current of a motor driving the pump. These features 
present in the stator frequency spectrum of the motor stator current may be caused by 

1 0 load effects of the pump on the motor rather than changes in the motor itself. The present 
invention provides a system and method for extracting (e.g., synthesizing) the fault data 
directly from the instantaneous motor current data. This data relates not only to pump 
machinery conditions, but also pump process conditions. Thus, by employing current 
signature analysis of the instantaneous current of the motor driving the pump, problems 

1 5 with the pump and/or process line can be detected without using invasive and expensive 
pressure and flow meters. Instead, a lower cost current sensor may be used and this 
sensor may be located in a motor control center or other suitable location remote from the 
motor and pump. 

Artificial neural networks (ANN) may thus be employed analyze the current 
20 signature data of the motor that relates to pump faults. Although, multi-iterative, 

supervised learning algorithms could be used, which could be trained and used only when 
a fully-labeled data set corresponding to all possible operating conditions, the application 
of unsupervised ANN techniques that can learn on-line (even in a single iteration) may be 
provided in accordance with the present invention. The current signature analysis, 
25 moreover, may be performed both on the pump to determine the operating state of the 
pump, and on the motor driving the pump so as to determine the operating state of the 
motor, simultaneously. 

The present invention also provides for preprocessing of the fault signature data 
before it is being used to train an ANN or design a decision module based on ANN 
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paradigms. The preprocessing eliminates outliers and performs scaling and bifurcation of 
the data into training and testing sets. Furthermore, it is desired to further post process 
the output generated by unsupervised ANN based decision modules for condition 
monitoring applications. This is because unsupervised ANN based decision modules 
5 when presented with a new operating condition can only signal the formation of a new 
output entry indicating that a possible new condition has occurred, but is not necessarily 
able to provide particular fault information. Post processing is carried out by utilizing the 
domain knowledge of a human expert to develop an expert system, or by correctly 
classifying this new operating state and encoding this information in a fuzzy expert 

1 0 system for future reference. 

Another aspect of the invention provides systems and methodologies for detecting 
motor faults by space vector angular fluctuation, which require no human intervention or 
downtime, in order to identify motor faults such as stator faults, rotor faults, and even 
imbalances in power applied to the motor in a timely fashion. Systems and 

15 methodologies are provided for detecting faults and adverse conditions associated with 
electric motors. The methodology provides for analyzing the angular fluctuation in a 
current space vector in order to detect one or more faults associated with the motor. 
Systems are disclosed having a diagnostics component adapted to obtain a space vector 
from a current signal relating to operation of the motor, and to analyze the space vector 

20 angular fluctuation in order to detect motor faults. 

To the accomplishment of the foregoing and related ends, the invention, then, 
comprises the features hereinafter fully described. The following description and the 
annexed drawings set forth in detail certain illustrative aspects of the invention. 
However, these aspects are indicative of but a few of the various ways in which the 

25 principles of the invention may be employed. Other aspects, advantages and novel 

features of the invention will become apparent from the following detailed description of 
the invention when considered in conjunction with the drawings. 
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Brief Description of the Drawings 

Fig. 1 is a schematic diagram illustrating an exemplary method in accordance 
with an aspect of the present invention; 

Fig. 2 is a side elevation view illustrating a motorized system having an 
5 exemplary diagnostics and control system in accordance with another aspect of the 
invention; 

Fig. 3 is a schematic diagram illustrating an exemplary diagnostics system with a 
neural network in accordance with another aspect of the invention; 

Fig. 4 is a schematic diagram illustrating another exemplary diagnostics system 
10 with an expert system in accordance with another aspect of the invention; 

Fig. 5 is a schematic diagram illustrating another exemplary diagnostics system 
with a data fusion component in accordance with another aspect of the invention; 

Fig. 6 is a schematic diagram illustrating another exemplary diagnostics system 
with a data fusion component, and neural network, and an expert system in accordance 
1 5 with another aspect of the invention; 

Fig. 7 is schematic flow diagram illustrating a portion of an exemplary 
diagnostics system in accordance with the invention; 

Fig. 8 is a schematic diagram further illustrating the exemplary diagnostics system 
of Fig. 7; 

20 Fig. 9 is a schematic diagram illustrating an exemplary cavitation classification in 

accordance with the invention; 

Fig. 10 is a perspective schematic diagram illustrating an exemplary neural 
network in accordance with another aspect of the invention; 

Fig. 1 1 is a function block diagram illustrating the diagnostics and control system 
25 in accordance with the present invention; 

Fig. 12 is a function block diagram further illustrating the diagnostics and control 
system of the present invention; 

Fig. 13 is a flow diagram for performing a Fast Fourier Transformation on a 
conditioned signal received from a motor in accordance with the present invention; 
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Fig. 14a is a graph of a comparison of a Fast Fourier Transform signal 
representative of a normal condition of a centrifugal pump and a cavitation condition of 
the same pump in accordance with the present invention; 

Fig. 14b is a magnified window of the graph of Fig. 14a in accordance with the 
5 present invention; 

Fig. 14c is a graph of a comparison of a Fast Fourier Transform signal 
representative of a normal condition of a centrifugal pump and a leakage condition of the 
same pump in accordance with the present invention; 

Fig. 14d is a magnified window of the graph of Fig. 14c in accordance with the 
1 0 present invention; 

Fig. 14e is a graph of a comparison of a Fast Fourier Transform signal 
yp representative of normal condition of a centrifugal pump and a faulty impeller condition 

JJ? of the same pump in accordance with the present invention; 

PI Fig. 14f is a magnified window of the graph of Fig. 14e in accordance with the 

y. 15 present invention; 

^ Fig. 14g is a table diagram of current signal amplitudes over a range of 

O frequencies, which may be employed to facilitate diagnosing an operating state of a 

;S J machine in accordance with the present invention; 

Fig. 1 5 is a flow diagram for creating a fault signature from the current spectrum 
M* 20 on the conditioned signal received from a motor in accordance with the present invention; 

Fig. 16 is a block diagram of an artificial neural network in accordance with the 
present invention; 

Fig. 17 is a block diagram of an artificial neural network in accordance with the 
present invention utilizing an Adaptive Resonance Theory paradigm; 
25 Fig. 18 is a block diagram of an artificial neural network in accordance with the 

present invention utilizing the second version of the Adaptive Resonance Theory 
paradigm; 

Fig. 19 is a functional schematic diagram of the diagnostics and control system 
including a stand-alone decision module adapted for pump diagnosis in accordance with 
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the present invention; 

Fig. 20 is a flow diagram for performing an adaptive preprocessing act for the 
stand-alone decision module of Fig. 19 in accordance with the present invention; 

Fig. 21 is a schematic representation of the Fuzzy Rule Base Expert System in 
5 accordance with the present invention; 

Fig. 22 is a schematic representation of the stand-alone decision module in 
accordance with the present invention; 

Fig. 23 is an example of a set of fuzzy rules in accordance with the present 
invention; 

1 0 Fig. 24 is a schematic representation of a system for diagnosing and controlling a 

plurality of pumps in accordance with the present invention; 
m Fig. 25 is a function block diagram illustrating the diagnostics and control system 

X utilizing wavelets in accordance with the present invention; 

^ Fig. 26 is an exemplary plot illustrating sampled three phase currents with 

U 15 unbalanced harmonics in accordance with the invention; 

j Fig. 27 is an exemplary plot illustrating sampled current space vectors for a 
O balanced and an unbalanced system in accordance with the invention; 

flj Fig. 28 is an exemplary plot illustrating space vector angular fluctuations in the 

%1 time domain in accordance with the invention; 

20 Fig. 29 is an exemplary plot illustrating frequency spectrum of space vector 

angular fluctuation with fault indicative frequencies in accordance with the invention; 

Fig. 30 is an exemplary plot illustrating an exemplary 2fs real time component 
obtained via a Goertzel algorithm in accordance with the invention; 

Fig. 3 1 is and exemplary plot illustrating a fluctuation in amplitude of a 2fs 
25 component due to stator fault in accordance with the invention; 

Fig. 32 is an exemplary plot illustrating another exemplary 2fs real time 
component obtained via a Goertzel algorithm in accordance with the invention; 

Fig. 33 is an exemplary plot illustrating an amplitude fluctuation in 2fs and 4fs 
components due to rotor resistance imbalance according to the invention; 
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Fig. 34 is an exemplary plot illustrating an amplitude fluctuation in 2fs and 4fs 
components due to a broken rotor bar in accordance with the invention; 

Fig. 35 is an exemplary plot illustrating sidebands due to rotor resistance 
imbalance in accordance with the invention; 

Fig. 36 is an exemplary plot illustrating oscillation of 2fs and 2sfs components 
due to rotor asymmetry in accordance with the invention; 

Fig. 37 is an exemplary plot illustrating a current space vector spectrum in 
accordance with the invention; 

Fig. 38 is an exemplary plot illustrating a voltage space vector spectrum in 
accordance with the invention; 

Fig. 39 is an exemplary plot illustrating a space vector angular fluctuation 
spectrum with components due to rotor imbalance in accordance with the invention; 

Fig. 40 is an exemplary plot illustrating an amplitude fluctuation in (l-s)2fs and 
(l+s)2fs components due to a broken rotor bar in accordance with the invention; 

Fig. 41 is an exemplary plot illustrating oscillation of a 2fs component due to a 
broken rotor bar in accordance with the invention; 

Fig. 42 is an exemplary plot illustrating a current space vector spectrum in the 
presence of a broken rotor bar in accordance with the invention; 

Fig. 43 is an exemplary plot illustrating a voltage space vector spectrum in the 
presence of a broken rotor bar in accordance with the invention; and 

Fig. 44 is an exemplary plot illustrating a space vector angular fluctuation 
spectrum in the presence of a broken rotor bar in accordance with the invention. 

Detailed Description of the Invention 

The various aspects of the present invention will now be described with reference 
to the drawings, wherein like reference numerals are used to refer to like elements 
throughout. The invention provides a diagnostics and control system for controlling a 
motorized system and diagnosing the health thereof, with a controller operatively 
associated with the motorized system and adapted to operate the motorized system in a 
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controlled fashion, and a diagnostics system operatively associated with the motorized 
system and adapted to diagnose the health of the motorized system according to a 
measured attribute associated with the motorized system. 

Referring initially to Fig. 1, an exemplary method 2 is illustrated for controlling a 
5 motorized system and diagnosing the health thereof. The method 2 comprises operating a 
motor in the motorized system in a controlled fashion, and diagnosing the health of the 
motorized system according to a measured attribute associated with the motorized 
system. While the exemplary method 2 is illustrated and described herein as a series of 
blocks representative of various events and/or acts, the present invention is not limited by 

10 the illustrated ordering of such blocks. For instance, some acts or events may occur in 
different orders and/or concurrently with other acts or events, apart from the ordering 
illustrated herein, in accordance with the invention. In addition, not all illustrated blocks, 
events, or acts, may be required to implement a methodology in accordance with the 
present invention. Moreover, it will be appreciated that the exemplary method 2 and 

1 5 other methods according to the invention, may be implemented in association with the 
motorized systems illustrated and described herein, as well as in association with other 
systems and apparatus not illustrated or described. 

Beginning at 4, the exemplary method 2 comprises measuring an attribute at 5, 
wherein the attribute is associated with a motorized system (e.g., motorized pump, fan, 

20 conveyor system, compressor, gear box, motion control device, screw pump, and mixer, 
hydraulic or pneumatic machine, or the like). The attribute measured at 5 may comprise, 
for example, vibration, pressure, current, speed, and/or temperature associated with the 
motorized system. At 6, the health of the motorized system is diagnosed according to the 
measured attribute. A diagnostics signal is provided at 7, which may be indicative of the 

25 diagnosed motorized system health, whereby the motorized system is operated at 8 

according to a setpoint and/or the diagnostics signal generated at 7. The provision of the 
diagnostics signal at 7 may comprise obtaining a frequency spectrum of the measured 
attribute and analyzing the frequency spectrum in order to detect faults, component wear 
or degradation, or other adverse condition in the motorized system, whether actual or 
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anticipated. The diagnosis may further comprise analyzing the amplitude of a first 
spectral component of the frequency spectrum at a first frequency. 

In order to provide the diagnostics signal, moreover, the invention may provide 
the measured attribute to a neural network, an expert system, a fuzzy logic system, and/or 
5 a data fusion component, which generates the diagnostics signal indicative of the health 
of the motorized system. For example, such frequency spectral analysis may be 
employed at 6 in order to determine faults or adverse conditions associated with the 
system or components therein {e.g., motor faults, unbalanced power source conditions, 
etc.). In addition, the diagnosis at 6 may identify adverse process conditions, such as 
10 cavitation in a motorized pumping system. The diagnostics signal generated at 7 may be 
^ employed by a controller associated with the system, whereby modified operation thereof 

i| may be performed in order to ameliorate or avoid actual or anticipated diagnosed health 

pi problems. 

% Another aspect of the invention provides systems and apparatus for controlling 

H 1 5 and diagnosing the health of motorized systems. Various aspects of the invention will be 
J" hereinafter illustrated with respect to an exemplary motorized pumping system. 

M However, it will be appreciated by those skilled in the art that the invention finds 

III application in association with motorized systems in addition to those illustrated and 

5| described herein, including but not limited to motorized pumps, fans, conveyor systems, 

^ 20 compressors, gear boxes, motion control devices, screw pumps, mixers, hydraulic or 
pneumatic machines, or the like. 

Referring now to Figs. 2-6, a motorized pump system 12 is illustrated having an 
exemplary diagnostics and control system 66 for controlling the system 12 and 
diagnosing the health thereof, in accordance with the present invention. As illustrated 
25 and described in greater detail hereinafter, the diagnostics and control system 66 
comprises a controller 71 operatively associated with the motorized system 12 and 
adapted to operate the system 12 in a controlled fashion, and a diagnostics system 70 
operatively associated with the motorized system 12 and adapted to diagnose the system 
health according to a measured attribute. The diagnostics and control system 66 of the 
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invention may measure attributes such as vibration, pressure, current, speed, and/or 
temperature in order to operate and diagnose the health of the system 12. 

An exemplary motorized pumping system 12 is illustrated in Fig. 2 having a 
pump 14, a three-phase electric motor 16, and a control system 18 for operating the 
5 system 12 in accordance with a setpoint 19. Although the exemplary motor 16 is 
illustrated and described herein as a polyphase synchronous electric motor, the various 
aspects of the present invention may be employed in association with single-phase motors 
as well as with DC and other types of motors. In addition, while the exemplary pump 14 
may comprise a centrifugal type pump, the invention finds application in association with 
1 0 other pump types not illustrated herein, for example, positive displacement pumps. The 
n control system 18 operates the pump 14 via the motor 16 according to the setpoint 19 and 

*? one or more measured process variables, in order to maintain operation of the system 12 

CP commensurate with the setpoint 19 and within the allowable process operating ranges 

j* specified in setup information 68. For example, it may be desired to provide a constant 

^ 15 fluid flow, wherein the value of the setpoint 19 is a desired flow rate in gallons per 
s minute (GPM) or other engineering units. 

|3 The pump 14 comprises an inlet opening 20 through which fluid is provided to the 

5 if pump 14 in the direction of arrow 22 as well as a suction pressure sensor 24, which 

O senses the inlet or suction pressure at the inlet 20 and provides a corresponding suction 

^ 20 pressure signal to the control system 1 8. Fluid is provided from the inlet 20 to an 

impeller housing 26 including an impeller (not shown), which rotates together with a 
rotary pump shaft coupled to the motor 16 via a coupling 28. The impeller housing 26 
and the motor 16 are mounted in a fixed relationship with respect to one another via a 
pump mount 30, and motor mounts 32. The impeller with appropriate fin geometry 
25 rotates within the housing 26 so as to create a pressure differential between the inlet 20 
and an outlet 34 of the pump. This causes fluid from the inlet 20 to flow out of the pump 
14 via the outlet or discharge tube 34 in the direction of arrow 36. The flow rate of fluid 
through the outlet 34 is measured by a flow sensor 38, which provides a flow rate signal 
to the control system 18. 
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In addition, the discharge or outlet pressure is measured by a pressure sensor 40, 
which is operatively associated with the outlet 34 and provides a discharge pressure 
signal to the control system 18. It will be noted at this point that although one or more 
sensors (e.g., suction pressure sensor 24, discharge pressure sensor 40, outlet flow sensor 
5 38, and others) are illustrated in the exemplary system 12 as being associated with and/or 
proximate to the pump 14, that such sensors may be located remote from the pump 14, 
and may be associated with other components in a process or system (not shown) in 
which the pump system 12 is employed. Alternatively, flow may be approximated rather 
than measured by utilizing pressure differential information, pump speed, fluid 
10 properties, and pump geometry information or a pump model. 
™ Alternatively or in combination, inlet and/or discharge pressure values may be 

II estimated according to other sensor signals and pump / process information. The system 
m 12 further includes an atmospheric pressure sensor 31, a pump temperature sensor 33 

(e.g., thermocouple, RTD, etc.), and a vibration sensor 37 (e.g., accelerometer or the 
iU 15 like), providing atmospheric pressure, pump temperature, and pump vibration signals, 
; respectively, to the control system 18. The invention finds application in association with 

!y motorized systems having fewer, more, or different combinations of sensors, apart from 

III the sensors illustrated and described herein, wherein sensors providing signals indicative 
S4 of other system variables may be employed in order to diagnose and control a motorized 
^ 20 system in accordance with the present invention. 

In addition, it will be appreciated that while the motor drive 60 is illustrated in the 
control system 18 as separate from the motor 16 and from an exemplary diagnostics and 
control system 66, that some or all of these components may be integrated. Thus, for 
example, an integrated, intelligent motor may include the motor 16, the motor drive 60 
25 and the diagnostics and control system 66. Furthermore, the motor 16 and the pump 14 
may be integrated into a single unit (e.g., having a common shaft wherein no coupling 28 
is required), with or without integral control system (e.g., control system 18, comprising 
the motor drive 60 and the diagnostics and control system 66) in accordance with the 
invention. 
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The control system 18 further receives process variable measurement signals 
relating to motor (pump) rotational speed, motor temperature, and motor vibration via a 
speed sensor 46, a motor temperature sensor 47, and a motor vibration sensor 48, 
respectively. As illustrated and described further hereinafter, a diagnostics system 70 
within the diagnostics and control system 66 may advantageously detect and/or diagnose 
actual or anticipated wear, degradation, failure, and/or faults associated with components 
of the motorized system 12 (e.g., and/or of components in a larger system of which the 
motorized system 12 is a part) as well as other system performance conditions, such as 
cavitation, blockage, or the like. 

For instance, the diagnostics system 70 may advantageously be employed in order 
to detect and/or diagnose cavitation in the pump 14 using a neural network classifier 
receiving suction and discharge pressure signals from sensors 24 and 40, respectively, as 
well as flow and pump speed signals from the flow and speed sensors 38 and 46. The 
motor 16 provides rotation of the impeller of the pump 14 according to three-phase 
alternating current (AC) electrical power provided from the control system via power 
cables 50 and a junction box 52 on the housing of the motor 16. The power to the pump 
14 may be determined by measuring the current provided to the motor 16 via a current 
sensor 49 and computing pump power based on current, speed, and motor model 
information. This may alternatively or in combination be measured and computed by a 
power sensor (not shown), which provides a signal related thereto to the control system 
1 8. Furthermore, the motor drive 60 may provide motor current, voltage, and/or torque 
information to the diagnostics and control system 66, for example, wherein pump input 
power information may be calculated according to the torque and possibly speed 
information. 

The control system 18 also comprises a motor drive 60 providing three-phase 
electric power from an AC power source 62 to the motor 16 via the cables 50 in a 
controlled fashion (e.g., at a controlled frequency and amplitude) in accordance with a 
control signal 64 from the diagnostics and control system 66. The diagnostics and 
control system 66 receives the process variable measurement signals from the suction 
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pressure sensor 24, the discharge pressure sensor 40, the flow sensor 38, and the speed 
sensor 46, together with the setpoint 19 and/or other sensor signals, and provides the 
control signal 64 to the motor drive 60 in order to operate the pump system 12 
commensurate with the setpoint 19 and/or a diagnostics signal 72 from a diagnostics 
5 system 70. 

In this regard, the diagnostics and control system 66 may be adapted to control the 
system 12 to maintain a desired fluid flow rate, outlet pressure, motor (pump) speed, 
torque, suction pressure, or other performance characteristic. Setup information 68 may 
be provided to the diagnostics and control system 66, which may include operating limits 

10 (e.g., min/max speeds, min/max flows, min/max pump power levels, min/max pressures 
allowed, NPSHR values, and the like), such as are appropriate for a given pump 14, 
motor 16, and piping and process conditions. 

The diagnostics and control system 66 comprises a diagnostics system 70, which 
is adapted to detect and/or diagnose cavitation in the pump 14, according to an aspect of 

15 the invention. Furthermore, the diagnostics and control system 66 selectively provides 
the control signal 64 to the motor drive 60 via a controller component 71 (e.g., which 
may implement one or more control strategies, such as proportional, integral, derivative 
(PID) control, or the like) according to the setpoint 19 (e.g., in order to maintain or 
regulate a desired flow rate) and/or a diagnostics signal 72 from the diagnostics system 

20 70 according to detected cavitation in the pump, whereby operation of the pumping 
system 12 may be changed or modified according to the diagnostics signal 72. The 
diagnostics system 70 as well as the controller 71 may be implemented in hardware, 
software, and/or combinations thereof according to appropriate coding techniques in 
order to implement the various aspects of the present invention. 

25 The diagnostics system 70 may detect the existence of one or more actual or 

anticipated conditions associated with the pump system 12, such as system component 
wear, degradation, failures, faults, or the like. In addition, the diagnostics system 70 may 
detect and diagnose process conditions associated with the system 12. For instance, the 
diagnostics system 70 may identify or detect cavitation in the pump 14, and additionally 
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diagnose the extent of such cavitation according to pressure and flow data from the 
sensors 24, 40, and 38 (e.g., and pump speed data from the sensor 46), or alternatively 
from current information from sensor 49, whereby the diagnostics signal 72 is indicative 
of the existence and extent of cavitation in pump 14. 

Referring also to Figs. 3-6, the diagnostics system 70 may comprise a neural 
network 80, an expert system 82, a data fusion component 84, and/or a fast Fourier 
transform (FFT) system 86 for generating a frequency spectrum 88, in order to provide 
the diagnostics signal 72 according to one or more measured system attributes (e.g., via 
signals from the sensors 24, 33, 35, 37, 38, 40, 46, 47, 48, and/or 49). The measured 
attributes may also be provided to the controller 71 from one or more of the sensors 24, 
33, 35, 37, 38, 40, 46, 47, 48, and/or 49 for performing closed loop control of the system 
12 in accordance with the setpoint 19 and/or the diagnostics signal 72. Although the 
aspects of the invention are described hereinafter within the context of employing a 
neural network 80, it is to be appreciated the invention may employ other nonlinear 
training systems/methodologies, for example, back-propagation, Baysian, Fuzzy Set, Non 
Linear regression, or other neural network paradigms including mixture of experts, 
cerebellar model arithmetic computer (CMACS), Radial Basis Functions, directed search 
networks, and functional link nets. 

The diagnostics and control system 66 thus comprises the controller 71 
operatively associated with the motorized system 12 to operate the system 12 in a 
controlled fashion via the control signal 64 to the motor drive 60, as well as the 
diagnostics system 70 operatively associated with the motorized system 12 and adapted 
to diagnose the health thereof according to one or more measured attributes (e.g., 
vibration, pressure, current, speed, and/or temperature), such as the signals from sensors 
24, 33, 35, 37, 38, 40, 46, 47, 48, and/or 49 in accordance with the present invention. It 
will be noted at this point that although the exemplary motorized system 12 includes a 
pump 14, that the diagnostics and control system 66 may be employed in association with 
other motorized systems (not shown) having a motor and a load, including but not limited 
to valve, pump, conveyor roller, fan, compressor, and/or gearbox type loads. 
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In accordance with the invention, the diagnostics system 70 provides diagnostics 
signal 72 according to the health of the system 12 and/or components therein (e.g., motor 
16, drive 60, pump 14), wherein the controller 71 may advantageously provide control 
signal 64 to the drive 60 according to the setpoint 19 and/or the diagnostics signal 72. As 
5 another application, the measured attribute or attributes may comprise vibration signals 
obtained from one or both of the sensors 37 and 48, and/or from other vibration sensors 
(e.g., accelerometers) associated with the motor 16 or the pump 14. For instance, the 
diagnostics system 70 may obtain one or more motor vibration signals (e.g., from sensor 
48), and diagnose the health of one or more motor bearings, motor shaft alignment (e.g., 

10 axial and/or angular alignment of the motor 16 with the pump 14 and/or misalignment 
thereof), and/or of the motor mounting (e.g. , via motor mounts 32) according to the 
measured vibration. 

The diagnostics system 70 in this regard, may be adapted to perform frequency 
spectral analysis of the measured vibration signal from sensor 48. Such spectral analysis 

15 may be performed, for example, via one or more of the neural network 80 and the expert 
system 82, where the diagnostics signal 72 indicates the health of the motorized system 
12 according to frequency spectral analysis of the measured vibration signal. 
Alternatively or in combination, the diagnostics system 70 may employ the data fusion 
system 84 in order to derive one or more vibration signals from at least one sensor 

20 associated with the motorized system 12. Thus, for example, where the vibration sensors 
37 and/or 48 are inoperative, data fusion techniques may be employed in accordance with 
the invention, to derive vibration information from other available system signals, such as 
via the current sensor 49. 

The present invention may thus employ data fusion in situations in order to take 

25 advantage of information fission which may be inherent to a process (e.g. , vibration in 
the motor 16) relating to sensing a physical environment through several different sensor 
modalities. In particular, one or more available sensing elements (e.g., sensors 24, 33, 
35, 37, 38, 40, 46, 47, 48, and/or 49) may provide a unique window into the physical 
environment where the phenomena to be observed is occurring (e.g., in the motorized 
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system 12 and/or in a system of which the motorized pumping system 12 is a part). 
Because the complete details of the phenomena being studied (e.g., detecting the 
operating state of the system 12 or components thereof) may not be contained within a 
single sensing element window, there is information fragmentation which results from 
this fission process. These information fragments associated with the various sensing 
devices may include both independent and dependent components. 

The independent components may be used to further fill out (or span) the 
information space and the dependent components may be employed in combination to 
improve the quality of common information recognizing that all sensor data may be 
subject to error and/or noise. In this context, data fusion techniques employed in the data 
fusion system 84 may include algorithmic processing of sensor data (e.g., from one or 
more of the sensors 24, 33, 35, 37, 38, 40, 46, 47, 48, and/or 49) in order to compensate 
for the inherent fragmentation of information because a particular phenomena may not be 
observed directly using a single sensing element. Thus, the data fusion system 84 
provides a suitable framework to facilitate condensing, combining, evaluating and 
interpreting the available sensed information in the context of the particular application. 
It will further be appreciated that the data fusion system 84 may be employed in the 
diagnostics and control system 66 in order to employ available sensors to infer or derive 
attribute information not directly measurable, or in the event of sensor failure. 

Thus, the present invention provides a data fusion framework and algorithms to 
facilitate condensing, combining, evaluating and interpreting various sensed data. The 
present invention also facilitates establishing a health state of a system employing the 
diagnostics and control system 66, as well as for predicting or anticipating a future state 
of the system 12 (e.g., and/or of a larger system of which the motorized pump system 12 
is a part). As illustrated in Fig. 6, the data fusion system 84 may be employed to derive 
system attribute information relating to any number of attributes according to measured 
attribute information (e.g., from the sensors 24, 33, 35, 37, 38, 40, 46, 47, 48, and/or 49) 
in accordance with the present invention. In this regard, the available attribute 
information may be employed by the data fusion system to derive attributes related to 
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failed sensors, and/or to other performance characteristics of the system 12 for which 
sensors are not available. Such attribute information derived via the data fusion system 
84 may be employed in generating the diagnostics signal 72, and/or in performing control 
functions in the controller 71. 
5 In another example, the measured attributes may comprise flow and pressure 

signals obtained from sensors 24, 38, and/or 40 associated with the pump 14, wherein the 
diagnostics system 70 provides the diagnostics signal 72 indicative of pump cavitation 
according to the measured flow and pressure signals. The invention thus provides for 
health indications relating to component conditions (e.g., wear, degradation, faults, 

10 failures, etc.), as well as those relating to process or systems conditions, such as 

cavitation in the pump 14. The diagnostics system 70 may comprise a classifier system, 
such as the neural network 80, detecting pump cavitation according to the measured flow 
and pressure signals, which may be provided as inputs to the neural network 80. The 
cavitation indication in the resulting diagnostics signal 72 may further be employed to 

15 modify operation of the system 12, for example, in order to reduce and/or avoid such 
cavitation. Thus, an appropriate control signal 64 may be provided by the controller 71 
to the motor drive 60 in order to avoid anticipated cavitation, based on the diagnostics 
signal 72 (e.g., and/or the setpoint 19), whereby the service lifetime of one or more 
system components (e.g., pump 14) may be extended. 

20 In another related example, cavitation (e.g., actual or suspected) in the pump 14 

may be detected via measured (e.g., or derived) current signal measurements, for 
example, via the sensor 49. The diagnostics system 70 in this instance may provide a 
diagnostics signal 72 indicative of pump cavitation according to the measured current. In 
order to detect cavitation using such current information, the diagnostics system 70 may 

25 employ the neural network 80 to synthesize a change in condition signal from the 
measured current. In addition, the diagnostics system 70 may further comprise a 
preprocessing portion (not shown) operatively coupled to the neural network 80, which 
conditions the measured current prior to inputting the current into the neural network 80, 
as well as a post processing portion operatively coupled to the neural network 80 to 
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determine whether the change in condition signal is due to a fault condition related to the 
motorized system 12. In this regard, the post processing portion may comprise a fuzzy 
rule based expert system, such as system 82. In addition, the diagnostics system 70 may 
detect one or more faults relating to the operation of the pump 14 and/or one or more 
faults relating to the operation of the motor 16 driving the pump 14 according to the 
measured current. 

Other faults may be detected and diagnosed using the diagnostics and control 
system 66 of the invention. For instance, the diagnostics system 70 may be adapted to 
obtain a space vector angular fluctuation from a current signal (e.g., from the current 
sensor 49) relating to operation of the motor 16, and further to analyze the space vector 
angular fluctuation in order to detect at least one fault in the motorized system 12. Such 
faults may include, for example, stator faults, rotor faults, and/or an imbalance condition 
in the power applied to the motor 16 in the motorized system 12 

In this situation, the diagnostics system 70 may obtain a current signal associated 
with the motor 16 from sensor 49, and calculate a space vector from the current signal. 
The diagnostics system 70 determines a space vector angular fluctuation from the space 
vector, and analyzes the space vector angular fluctuation in order to detect one or more 
faults associated with the motor 16. For instance, first, second, and third phase current 
signals associated with the motorized system 12 may be sampled in order to obtain the 
current signal, and corresponding first, second, and third phase space vectors may be 
computed in the diagnostics system 70. 

A resulting space vector may then be calculated, for example, by summing the 
first, second, and third phase space vectors. The diagnostics system 70 may then 
compare the space vector with a reference space vector, wherein the reference space 
vector is a function of a constant frequency and amplitude, and compute angular 
fluctuations in the space vector according to the comparison, in order to determine the 
space vector angular fluctuation. The diagnostics system 70 then performs frequency 
spectrum analysis {e.g., using the FFT component 86) of the space vector angular 
fluctuation to detect faults associated with the motorized system 12. For example, motor 
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faults such as rotor faults, stator faults, and/or unbalanced supply power associated with 
the pump motor 16 may be ascertained by analyzing the amplitude of a first spectral 
component of the frequency spectrum at a first frequency, wherein the diagnostics system 
70 may detect fluctuations in amplitude of the first spectral component in order to detect 
one or more faults or other adverse conditions associated with the motorized system 12. 

In this regard, certain frequencies may comprise fault related information, such as 
where the first frequency is approximately twice the frequency of power applied to the 
motor 16. Alternative to generating a full spectrum, the diagnostics system 70 may 
advantageously employ a Goertzel algorithm to extract the amplitude of the first spectral 
component in order to analyze the amplitude of the first spectral component. The 
diagnostics signal 72 indicating such motor faults may then be employed by the 
controller 71 to modify operation of the pumping system 12 to reduce or mitigate such 
faults. 

Referring now to Fig. 7, the processing performed on vibration and other sampled 
data from one or more of the system sensors 24, 33, 35, 37, 38, 40, 46, 47, 48, and/or 49 
in the diagnostics system 70 may comprise demodulation. One demodulation technique, 
sometimes referred to as enveloping, may be performed by the system 70 in order to 
synthesize sampled digital attribute data 100 into a form usable for detecting and 
diagnosing component degradation and process conditions in the system 12. The 
demodulation is illustrated in Fig. 7 in association with vibration data (e.g., from one of 
the sensors 37 or 48). The digital vibration data 100 enters the diagnostics system 70 and 
passes through a band pass filter 102, which removes frequencies outside the scope of 
interest and within the dynamic range of the system 70 to form a filtered signal 104. 

The filtered signal 104 passes through a rectifier 106, for example a diode, which 
forms a rectified signal 108. The rectified signal 108 passes through a low pass filter 1 10 
which removes the high frequencies to form a relatively low frequency signal 1 12. The 
low frequency signal 1 12 is passed through a capacitor 1 14 to produce a demodulated 
signal 1 1 6. A fast Fourier transform (FFT) is performed on the demodulated signal 1 1 6 
by FFT operator 118 (e.g., such as FFT component 86 of Fig. 6) in order to produce a 
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vibration spectrum 120, for example, using commercially available fast Fourier transform 
software such as MATLAB by The Math Works. The FFTs of the vibration signal data 
are discretized over N number of points to facilitate processing, such as where N=2,048, 
however, it will be appreciated that the FFTs of each signal may be discretized over any 
5 suitable number of such points. The vibration spectrum 1 20 may then be analyzed by a 
host processor in the control system 18 in order to determine the health of the motor 16 or 
other component in the motorized system 12. 

Although the demodulation has been described with respect to obtaining FFTs of 
vibration or other signals, other suitable techniques may be employed. For example, 
1 0 wavelet transforms may be taken of the sensor data. One advantage to using the wavelet 
transform is that the total size of the transform is a compact representation of the original 
signal and will require considerably less storage space than the original signal. 
Alternatively, where particular spectral components are of interest, techniques such as the 
Goertzel algorithm may be used to obtain such components without the need to generate 
1 5 an entire frequency spectrum. Furthermore, such spectral analysis of vibration data to 

detect and diagnose bearing faults, may further take into consideration motor (e.g. , pump) 
speed information (e.g., from sensor 46 in order to suitably interpret the FFT and detect 
defects and isolate particular bearing faults. 

Referring also to Fig. 8, the diagnostics system 70 may further comprise a pre- 
20 processing component 202 receiving pressure and flow data from the sensors 24, 40, and 
38, respectively, which provides one or more attributes 204 to the neural network 80, 
wherein the attributes 204 may represent information relevant to cavitation in the pump 
14. The attributes 204 may be extracted from the measured pressure, flow, and/or speed 
values associated with the pumping system 12, and used to characterize pump cavitation 
25 by the neural network 80. The neural network 80, in turn, generates a diagnostics signal 
72 which may comprise a cavitation classification 206 according to another aspect of the 
invention. The neural network classifier 80 thus evaluates data measured in the 
motorized pumping system 12 (e.g., represented by the attributes 204) and produces a 
diagnosis (e.g., diagnostics signal 72) assessing the presence and severity of cavitation in 
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the system 12. The neural network 80 in this regard, may employ one or more 
algorithms, such as a multi-layer perception (MLP) algorithm in assessing pump 
cavitation. 

As illustrated further in Fig. 9, the diagnostics signal 72 output by the classifier 

5 neural network 80 is indicative of both the existence and the extent of cavitation in the 
pumping system 12. For instance, the exemplary signal 72 comprises a classification 206 
of pump cavitation having one of a plurality of class values, such as 0, 1, 2, 3, and 4. In 
the exemplary classification 206 of Fig. 9, each of the class values is indicative of the 
extent of cavitation in the pumping system 12, wherein class 0 indicates that no cavitation 

10 exists in the pumping system 12. The invention thus provides for detection of the 
existence of cavitation (e.g., via the indication of class values of 1 through 4 in the 
diagnostics signal 72), as well as for diagnosis of the extent of such detected cavitation, 
via the employment of the neural network classifier 80 in the diagnostics system 70. 
Referring now to Fig. 10, an exemplary neural network 80 comprises an input 

15 layer 210 having neurons 212, 214, 216, and 218 corresponding to the suction pressure, 
discharge pressure, flow rate, and pump speed signals, respectively, received from the 
sensors 24, 40, 38, and 46 of the pumping system 12. One or more intermediate or 
hidden layers 220 are provided in the network 80, wherein any number of hidden layer 
neurons 222 may be provided therein. The neural network 80 further comprises an 

20 output layer 230 having a plurality of output neurons corresponding to the exemplary 

cavitation classification values of the class 206 illustrated and described hereinabove with 
respect to Fig. 9. Thus, for example, the output layer 230 may comprise output neurons 
232, 234, 236, 238, and 240 corresponding to the class values 0, 1, 2, 3, and 4, 
respectively, whereby the neural network 80 may output a diagnostics signal (e.g., signal 

25 72) indicative of the existence as well as the extent of cavitation in the pumping system 
(e.g., system 12) with which it is associated. 

In this regard, the number, type, and configuration of the neurons in the hidden 
layer(s) 220 may be determined according to design principles known in the art for 
establishing neural networks. For instance, the number of neurons in the input and output 
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layers 210 and 230, respectively, may be selected according to the number of attributes 
{e.g., pressures, flow, speed, etc.) associated with the system 70, and the number of 
cavitation classes 206. In addition, the number of layers, the number of component 
neurons thereof, the types of connections among neurons for different layers as well as 
5 among neurons within a layer, the manner in which neurons in the network 80 receive 
inputs and produce outputs, as well as the connection strengths between neurons may be 
determined according to a given application {e.g. , motorized system) or according to 
other design considerations. 

Accordingly, the invention contemplates neural networks having many 
1 0 hierarchical structures including those illustrated with respect to the exemplary network 
80 of Fig. 10, as well as others not illustrated, such as resonance structures. In addition, 
the inter-layer connections of the network 80 may comprise fully connected, partially 
connected, feed-forward, bi-directional, recurrent, and off-center or off surround 
interconnections. The exemplary neural network 80, moreover, may be trained according 
15 to a variety of techniques, including but not limited to unsupervised learning, 

reinforcement learning, and back propagation, wherein the learning may be performed 
on-line and/or off-line. Furthermore, the training of the network 80 may be accomplished 
according to any appropriate training laws or rules, including but not limited to Hebb's 
Rule, Hopfield Law, Delta Rule, Kohonen's Learning Law, and/or the like, in accordance 
20 with the present invention. In addition, although one or more aspects of the present 

invention are primarily described herein in the context of employing a neural network, it 
is to be appreciated the invention may employ other nonlinear training systems and/or 
methodologies (e.g., example, back-propagation, Baysian, Fuzzy Set, Non Linear 
regression, or other neural network paradigms including mixture of experts, cerebellar 
25 model arithmetic computer (CMACS), Radial Basis Functions, directed search networks, 
functional link nets, and the like). 

Referring now to Fig. 1 1, the diagnostic system 70 may further perform one or 
more types of signal conditioning on the measured attribute signals from the system 
sensors 24, 33, 35, 37, 38, 40, 46, 47, 48, and/or 49. As an example, one or more motor 
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currents may be sampled, for example, from the sensor 49, and processed accordingly, in 
order to detect and/or diagnose cavitation conditions in the pump 14. In this regard, the 
current sensor 49 may comprise a current transformer providing a current signal to a 
signal conditioning system 306, which in turn provides a conditioned current signal 320 
5 to an analog to digital (A/D) converter 272. Thereafter, a digitized current signal is 
provided from the A/D 272 to a processor 290 in the diagnostics and control system 66 
for further processing, as illustrated in Fig. 1 1 . 

Fig. 1 1 illustrates the functional events or acts that the current data from the motor 
16 driving the pump 14 is subjected to for pump condition diagnostics. After the current 
1 0 data from current transformer 49 is conditioned by signal conditioning circuit 3 06, it is 
converted from analog data to digital data by the A/D converter 272, so that it can be 
further processed by processor 290. Processor 290 first performs the acts of computing 
Fast Fourier Transforms 326 of the current data. The processor 290 controls the signal 
sampling and digitizing rate as well as any buffering of the digitized signals that might be 
1 5 needed. This data collection rate should be selected to provide sufficient data upon 

which the processor 290 can generate a comprehensive frequency spectrum of the motor 
current signal suitable for analysis using commercially available Fast Fourier Transform 
software, such as for example MATLAB by The Math Works. 

The spectral analysis basically involves the application of a 'Hanning Window' to 
20 the time-domain data, calculation of power spectrum from each set of the windowed 

time-domain data for the specified number of sets, and then finding the average spectrum 
using the Welch method. A flow chart of the general scheme is shown in Fig. 13. The 
output spectra corresponding to the A/D current signal, is stored into a designated output 
file for future use. The parameters for the A/D data acquisition, such as the number of 
25 sets (NumSets), number of samples per set (NumPts) and the sampling rate are selected 
to be 8, 8192 and 4096 Samples/sec respectively, however it will be appreciated that any 
appropriate values for these parameters may be employed in accordance with the present 
invention. These parameters yield a frequency resolution of 0.5 Hz and a bandwidth of 0 
to 2048 Hz in the frequency spectra. The time-domain data consists of a contiguous 
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record of 65536 data values collected over 16 seconds, which is then divided into eight 
equal sets. The noise smoothing is of a satisfactory level after the averaging of eight 
consecutive spectra. A Harming window may be employed for the windowing purpose 
because of its ability to reduce the leakage effect to a minimum. 
5 The processor 290 then finds and stores fault signatures relating to the operation 

of the pump in 328. Figs. 14a-14f show how pump faults can be detected comparing 
frequency spectrums of current motor data relating to normal and fault conditions. Figs. 
14a- 14b show both a frequency spectrum of a centrifugal pump in a normal condition 350 
and a cavitation condition 352. Fig. 14b illustrates a blown up version of Fig. 14a along a 

1 0 limited frequency range, wherein a difference in the spectrums can be seen during 
cavitation. Furthermore, Figs. 14c and 14d show both a frequency spectrum of a 
centrifugal pump in a normal condition 354 and a leakage condition 356, while Figs. 14e 
and 14f illustrate both a frequency spectrum of a centrifugal pump in a normal condition 
358 and a faulty impeller condition 360. It should be appreciated that other fault 

15 conditions such as blockage will also lead to different frequency spectrum characteristics 
of the current data of the motor driving the pump. Furthermore, different types of pumps 
will exhibit different types of faults with associated different spectrum characteristics. 
For example, pump faults could include cavitation, blockage, two-phase flow, impeller 
wear, impeller damage, the impeller impacting with the casing, a pump out of balance, 

20 corrosion, surge/hammer or other types of faults. 

In one aspect of the invention, the processor 290 could access a table, such as the 
one shown in Fig. 14g. The table 400 is shown which the processor 290 accesses when 
performing signature analysis to diagnose the health of the pump 14. The table 400 
includes current amplitude data (Ao thru Az) over a range of frequencies (fo thru fz). The 

25 table 400 is stored in memory in the diagnostics and control system 66, so as to be easily 
accessible by the processor 290. The table 400 includes various health states of the pump 
shown generally at 402, which correspond to current amplitudes over the frequency range 
fo thru fn. For example, referring to the row identified by reference numeral 404, when 
the current amplitudes are A 2 34 at f 0 , A 27 at fi, A478 at f 2 , A24 at f 3 , A127 at f 4 , Aq at f n , 
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the table 400 indicates that the pump 14 has a "pump fault 6". The "pump fault 6" could 
be a cavitation fault or a variety of other pump related faults. As will be appreciated, the 
table 400 can store N number of current signatures corresponding to various health states 
of the pump 14, which the processor 290 can employ to diagnose the health of the pump 
5 14. 

Furthermore, certain discriminating fault attributes may be extracted from the 
frequency spectrum of motor current, which relate certain fault conditions of the pump. 
Typical attributes such as motor slip and noise directly relate to degree of cavitation. An 
algorithm performed by processor 290 can be seen in Fig. 1 5, which evaluates the 

10 components of the frequency spectrum and derives certain attributes for a centrifugal 
pump. The attributes shown in Fig. 15 are slip, FsAmp, SigAmp, Noise J, Noise_2, 
Noise_3, Noise_4 and Noise_5, which form the fault signature for the centrifugal pump. 
For example, if FsAmp is low and slip is low there is fault condition of severe blockage. 
Other types of faults can be found, such as cavitation and faulty impeller by evaluating 

15 the above attributes. The attributes then may be subjected to the act of preprocessing 330 
to be acceptable by the neural network 80. 

In some cases may be desirable to preprocess the fault signature data before it is 
being used to train or design a decision module based on ANN paradigms. The 
preprocessing can be divided into three tasks, namely, elimination of outliers, scaling, 

20 and bifurcation of data into training and testing sets. Elimination of outliers is concerned 
with detecting any such data pattern that has one or more attributes which seem to have 
an extraordinarily large or small values compared to the allowed range for that 
attribute(s). Such data patterns are known as outliers, which could be generated due to 
errors during data collection or fault signature formation or due to noise. Elimination of 

25 such data patterns facilitates proper utilization of the data set for designing the network. 
The adverse effects caused by not eliminating outliers are compression of the useful 
range of a given attribute during scaling and causing difficulties for the network in 
converging to a final solution. 

The fault signature, which is an array of real values, is known as an analog data 
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pattern, analog exemplar, or feature vector in the field of neural networks. Although it is 
possible to take the attributes generated in the previous section and apply them directly as 
inputs to a neural network, it is practically more beneficial in terms of simplicity of the 
designed network and in terms of the computational load on the processor, to scale the 

5 attributes in such a way that each of the attributes has similar boundaries such as {0, 1 } 
or {-1, 1}. For example, in the pump condition monitoring data patterns, the slip attribute 
would have a value on the order of 10" 2 while the FsAmp attribute can have values 
greater than 10 3 and the possible range of their values can also differ greatly. Using such 
a data set directly, without scaling, in conjunction with a neural network would lead to 

10 large values of network parameters (e.g., weight vectors) and the time taken for 

completely training the network or designing the network would be high. The simplest 
and most widely used scaling method is the linear scaling method. Alternatively 
sigmoidal scaling can be used. 

Returning to Fig. 1 1, the processor 290 then transmits the fault signature to neural 

15 network 80 after 330. Neural network 80 could be a multi-iterative, supervised learning 
algorithm such as feed forward, back propagation network. This network could be 
trained and used with a fully labeled data set corresponding to all possible operating 
conditions. However, the application of unsupervised ANN techniques that could learn 
on-line (even in a single iteration) may be desirable. Also, the ability to perform 

20 incremental learning such as provided in the proposed unsupervised neural network is 

desirable. Processor 290 can be programmed to perform the necessary post processing of 
the fault data in 332 and provide that information in the diagnostics signal 72. It should 
be appreciated that the functions of diagnosing the motor current data containing pump 
fault information can be completely performed by the diagnostics and control system 66. 

25 Fig. 12 illustrates an implementation of the present invention wherein the 

diagnostics and control system 66 performs the signal conditioning 380, AID conversion 
382, spectral analysis 384, formation of the fault signatures 386 and pre-processing 388. 
The diagnostics and control system 66 may also include a decision module 390 which 
may be comprised of neural network 80 and post processing 332. 
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Neural network 80, which may comprise an Artificial Neural Network (ANN), 
will now be discussed with regards to another implementation in accordance with the 
present invention. ANNs are a set of algorithms or processing acts that can be used to 
impart the capabilities such as generalization and pattern recognition to a standard 
5 computer system. These algorithms can learn to recognize patterns or classify data sets 
by two methods, known as supervised learning and unsupervised learning. Considering 
the block diagram shown in the Fig. 16, in which an ANN is represented by a Square 
block' 410, the ANN is supervised if the network is provided with a set of input patterns 
412 along with their designated outputs 414, and the network learns by changing its 
1 0 internal parameters in such a way that it produces the corresponding designated output 
*»* pattern for a given input pattern. The ANN is known as unsupervised if the network 

V chooses its own output class for a given input pattern without any external supervision or 

:0i feed-back. 

% The ANN structure basically contains processing nodes called neurons arranged 

H; 15 in two or more layers. The nodes are extensively interconnected typically in a feed 
, s " forward manner through connections and associated interconnection weights. For a 

pattern recognition application, the input pattern could be a digitized image of an object 
tU to be recognized and the output could be the same image (e.g., with reduced noise), or it 

fl could be a class representation of the image. For a condition monitoring application such 

^ 20 as for motorized systems, the input pattern may include preprocessed fault signature data 
derived from the stator current spectra, and the output is a class representation of various 
fault conditions. There are various architectures and learning schemes for the ANNs and 
some of the ANNs may be more suitable for use as decision makers for condition 
diagnostics applications than others. 
25 It has been found that one-shot unsupervised ANN paradigms may be more 

suitable for the development of condition monitoring system than supervised ANN 
paradigms. Furthermore, a monitoring system based on an unsupervised neural network 
that can learn about the conditions of a process plant from a single pass of the training 
data, can provide a better solution for detecting new plant conditions. Supervised 
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networks must be trained off line, and thus have a fixed set of variables and cannot give a 
valid output upon receipt of a new condition. A block diagram of a one such one-shot 
unsupervised ANN network known as the Adaptive Resonance Theory (ART) is shown 
in Fig. 17. 

5 ART has two layers namely the input (or comparison) layer 440 and output (or 

recognition) layer 442, These layers are connected together, unlike the other networks 
discussed above, with feedforward (denoted by W) as well as feedback (denoted by T) 
connections. The neurons of the output layer also have mutual connections useful for 
lateral inhibition (not shown in the Figure). The signals Control- 1 and Control-2 are 

1 0 responsible for controlling the data flow through the input and output layers, 440 and 
442, respectively. The reset circuit is responsible for determining the effectiveness with 
which a winning output neuron represents the input pattern. It is also responsible for 
resetting the ineffective neurons and designating a new neuron for representing a given 
input pattern. The training of the ART network can be done in either a fast learning 

1 5 mode or in a slow learning mode. The fast learning mode allows the network's feed 

forward weights to attain their optimum values within few learning cycles (epochs) while 
the slow learning mode forces the weights to adapt over many epochs. The fast learning 
mode can be used to train the network even in a single epoch. This is appropriate for the 
present case because the salient features of the problem domain are already well defined 

20 in the fault signatures. 

Fig. 1 8 illustrates the ART-2 architecture, which is essentially the same as the 
ART architecture except that the input layer of ART-2 has six sublayers (w,x,v,u,p and q) 
which are designed to enhance the performance of the network to cope with the 
continuously varying real input data. The additional sublayers incorporate the effects of 

25 feature enhancement, noise suppression, sparse coding and expectation from the 
recognition layer. The reset function is incorporated by an additional sublayer (r) 
containing the same number of nodes as are in the sublayers of the input layer. 

Another useful unsupervised ANN algorithm applicable to the present invention is 
the Associative List Memory (ALM) paradigm. Associative List Memory (ALM) is an 
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autoassociative memory capable of handling long bit strings (up to 10 4 bits) of data. 
ALM's learning mechanism is based on unsupervised training technique and can learn 
within a single epoch. It has been successfully applied to reduce and classify the data 
produced by satellites and the results have shown that ALM is comparable or better than 
5 other associative memories such as the Sparse Distributed Memory for this application. 
It provides direct access to the learned exemplars and can be implemented on a low-cost 
computing platform (e.g., a 16 bit microprocessor). With suitable modifications ALM 
can be used for classifying analogue inputs {e.g., pump fault signature data). Its one-shot 
learning capability, simple algorithm structure, and incremental learning capability make 
10 it a suitable choice for implementing a stand-alone decision module for pump 
diagnostics. 

In a typical commercial or industrial motorized system, it may not be possible to 
generate a dataset containing all possible operating conditions or faults. Simple linear 
scaling of the available data may not ensure that any fixture data, which might occur due 

15 to a new plant operating condition, would lie within the same initial limits. Hence an 
adaptive preprocessing scheme may be employed in accordance with the invention, 
which adaptively updates the maximum and minimum values of each attribute so that the 
entire database is maintained to have values between 0 and 1. This increases the 
computational task, but the potential advantages of the generic condition monitoring 

20 scheme out- weigh this minor drawback. Output generated by an unsupervised neural 
network needs to be post processed in order to obtain useful classification information. 
In the classification applications using the unsupervised ANN paradigms, usually the 
domain expert (a human expert, having knowledge about the problem) examines the 
output classes assigned by the ANN in conjunction with system details such as which of 

25 the input patterns has activated which of the output entries, the class to which a given 
input actually belongs etc., and then assigns the class representation to each of the output 
entries. The above procedure is suitable for applications such as classification of satellite 
pictures etc., but may not be suitable for on-line condition monitoring applications. 
Without the above mentioned post processing, an unsupervised ANN based 
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decision module for the condition monitoring applications, can only signal the formation 
of a new output entry indicating that a possible new plant condition has occurred, but it 
cannot diagnose the new condition to be a particular case, such as a new normal 
condition, cavitation etc. It is, however, possible that an intelligent post processor based 
5 on expert system rules, can be used to assign the class representation to each of the output 
entries formed by the unsupervised neural networks such as ART and ALM. The 
attributes of the fault signatures obtained from the stator current of the motor are derived 
based on the physical interpretation of the functioning of the pumps and pumping 
systems. This means that the variation of each of the attributes are associated with the 
10 physical status of the pump. Thus, it is possible to correlate the attribute representing the 
pump condition with rules such as 'IF attribute- 1 is high and attribute-2 is high THEN the 
pump condition is cavitation', etc. An expert system based on such generic rules can be 
used to assign the class representation to the output entries formed by the unsupervised 
ANN paradigm so that the decision module produces a meaningful output without any 
1 5 external human supervision. 

Referring to Fig. 19, a condition monitoring system is provided that utilizes a 
stand-alone decision module 470 within the diagnostics and control system 66. The 
stand-alone decision module 470 includes hardware and/or software that performs an 
adaptive preprocessing act, an ALM neural network and software that utilizes a Fuzzy 
20 Expert System in a post processing act. To exploit the potential advantages of one-shot 
methods as decision modules for the on-line condition monitoring schemes, a flexible 
preprocessing method is employed which adaptively preprocesses the incoming data to 
bind its attribute values to the predefined limits of 0 to 1 , by using pre-specified 
maximum and minimum limits, which can be easily defined by the domain expert, who 
25 has knowledge of the equipment (type of motor, pump etc.) and the plant (flow rate, 
differential pressure head etc.). 

Where attribute values are not predictable, limits are set based on observation of 
the range encountered during measurement. These initial limits act as maximum and 
minimum boundaries for calculating scaling factors to preprocess the incoming data as 
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long as the incoming data patterns fall within range. When any of the attributes of the 
incoming data patterns has a value beyond these limits then the maximum and minimum 
boundaries are expanded to include the present pattern and new scaling factors are 
obtained. Whenever such a change of maxima or minima boundaries occurs, the neural 
5 network paradigm that uses the preprocessed data may be updated to account for the new 
changes in the scaling factors and also the previously accumulated data base may need to 
be rescaled using the new scaling factors based on how much the new boundaries have 
expanded. 

In the case of the ART-2 algorithm, the feed-forward and feed-back weight 

1 0 vectors may be rescaled to account for the new changes, and for an ALM the memory 
entries themselves rescaled. Although the proposed adaptive preprocessing scheme will 
work without defining any maximum and minimum limits beforehand, it is faster and 
more efficient to specify and utilize initial expected limits. It should be appreciated that 
the present invention could employ operator input training techniques. Fig. 20 illustrates 

15 a flow-chart of the functional description of the adaptive preprocessing scheme 472. The 
expected maxima and minima for each attribute are first derived by analyzing the normal 
condition data and the other relevant pump/plant details, ratings and physical system. It 
should be noted that in some circumstances it may not be necessary to employ 
preprocessing and postprocessing techniques. 

20 A block diagram of a Fuzzy Rule Based Expert System (FRBES) 476 is shown in 

Fig. 21 . It comprises three modules namely, a Fuzzy Rule Base 500, a Fuzzy Inference 
Engine 502 and a User Interface 504. The Fuzzy Rule Base 500 consists of Fuzzy type 
'IF - THEN' rules, while the Fuzzy Inference Engine 502 is the core logic program which 
includes details of the Fuzzy system such as the number and types of membership 

25 functions of inputs and outputs, their range, method of defuzzification, etc. The User 
Interface 504, in a production-level expert system is capable of querying the user for 
additional information in case the expert system is unable to reach a conclusion with the 
existing information in its data base. 

As can be seen in Fig. 22, the decision module 470 can be integrated as a stand- 
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alone decision module where the input can be received both by the Adaptive 
Preprocessor 472 and the FRBES 476. Bidirectional communication is set up between 
the ALM Neural Net 474 and the FRBES 476. The ALM may receive new classification 
information which may be entered from the user interface portion of the FRBES 476. 

5 The ALM Neural Net 474 then generates the appropriate output pattern to a user 

interface. The rule base of the FRBES in the present example comprises fifteen generic 
rules formed on the basis of the generic knowledge about the behavior of the attributes 
with respect to a given plant condition. These rules are listed in Fig. 23. TheMATLAB 
source code corresponding to the FRBES described above can be implemented using 

1 0 MATLAB's Fuzzy Logic Tool-Box. It should be appreciated that the Adaptive 

Preprocessor 472, ALM paradigm 474, and the FRBES 476 could be implemented in 
stand-alone hardware modules or in a single integrated hardware module within the 
system 66. Furthermore, separate processors and software programs could be used or a 
single processor and software program could be used to run the Adaptive preprocessing 

1 5 and FRBES acts. It should be appreciated that the stand-alone decision module 470 can 
be used to diagnose the operating condition of the pump system 12 (e.g., cavitation), the 
health or condition of the pump 14 (e.g., impeller damage), and also the condition of the 
motor 16 driving the pump 14 (e.g., rotor bar damage). 

Turning now to Fig. 24, the diagnostics and control system 66 is shown as part of 

20 a system 530 where the diagnostics and control system 66 performs classical current 
signature analysis on a plurality of motors 230 A -230 D . The motors 230 A -230 D drive 
pumps 232 A -232 D . Although only four motors and four pumps are shown as part of the 
system 530, it will be appreciated that virtually any number of motors and pumps may be 
employed in this system 530. The motors 230 A -230 D are each individually tied to the 

25 diagnostics and control system 66. The system 66 may include motor starters to start and 
stop the motors along with circuit breakers to protect the electric motors and electric 
wiring. Each respective motor 230 A -230 D has current sensor(s) coupled to its respective 
lead wire(s) to obtain current signal data. The digitized current signal data 320 for each 
motor 230 A -230 D is applied to a channel interface 532. The channel interface 532 
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includes a separate channel for each motor of the system 530. The channel interface 532 
is coupled to the diagnostics and control system 66. 

The diagnostics and control system 66 cycles through each channel at 
predetermined intervals to obtain current data with respect to each motor 230 A -230 D . The 
5 diagnostics and control system 66 then essentially carries out the same acts as described 
above with reference to Fig. 12 as to performing classical fault signature analysis to 
determine the operating condition of the pumps 232 A -232 D . The advantage of this system 
530 is that it allows for the analyzing a plurality of pumps by a single diagnostics and 
control system 66. In this way, a user of the diagnostics and control system 66 could 
1 0 monitor and analyze every motor or machine within a facility from a single location. 

Although the present invention has been described with respect to obtaining Fast 
Fourier Transforms of the current signals, it should be appreciated that other suitable 
techniques may be employed. For example, wavelet transforms may be derived from the 
current data. Fig. 25 shows replacing the act 326 (Fig. 1 1) of computing the Fast Fourier 
1 5 Transformation of digital conditioned signal 320 with computing the wavelet of the 
conditioned signal 320. All that is required for performing a wavelet transform is an 
appropriate set of analysis and synthesis filters. By using the wavelet transform with the 
neural network 80, a much smaller, compact, training set may be employed, which still 
enables the present invention to correctly classify the operating state of the motor pump 
20 system. Furthermore, original signal information from the wavelet coefficients as may be 
needed in the future. This approach also affords for a pseudo frequency domain and time 
domain analysis of the signal data. Such a combination involves relative simplicity of 
implementation while affording great flexibility in accommodating a broad range of 
signal types and noise levels. 
25 It should be appreciated that the present invention employs techniques that show 

direct application to machinery diagnosis using a variety of techniques, algorithms and 
models. These techniques could be readily expanded to include pump hardware and 
pump process protection via automatic shutdown, failure prediction/ prognostics, 
corrective action recommendations, monitor and control energy usage and to ensure EPA 
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and safety guidelines are complied with using monitoring and archival data storage. It 
should also be appreciated that the techniques provided in the present invention could be 
applicable to a broad range of pumps {e.g., centrifugal, positive displacement, 
compressors, vacuum pumps, etc.). 
5 In order to further illustrate the various aspects of the invention, Figs. 26-44 and 

the following description are provided, including exemplary experimental and simulated 
results showing the diagnostic capabilities according to the invention. The space vector 
angular fluctuation technique (SVAF), provides significant advantages over prior 
analytical methodologies, such as the zero cross times method (ZCT) for motor failure 

1 0 prediction. The ZCT method employs zero crossing times of three phase motor current 
waveforms as data for spectral analysis of induction motor current. The inventors of the 
present invention have found that fluctuations in the angle of the space vector hold 
information on motor condition, and that when these are analyzed by means of FFT or 
other frequency spectral analysis techniques, diagnostic indices for stator and rotor faults 

15 can be defined. The main drawback in the ZCT method was limited sampling frequency, 
whereas the SVAF method overcomes this limitation, thus giving more reliable 
diagnostic data, since aliasing effects are removed. 

The ZCT method of induction motor fault detection measures times t x at which the 
three phase currents cross through zero, A series of data values is derived as the time 

20 difference between pairs of adjacent zero crossing times minus the expected 60 degree 
time interval AT between two zero crossings. For a three phase system with six zero 
crossings per mains cycle, this produces six samples per supply cycle with the data given 
by 

25 gt^ti-tu-AT, (1) 

with AT = 1/300 sec for a 50 Hz supply. 

Since the phase lag angle of motor current behind the supply voltage varies with 
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load, it follows that fluctuations in load or speed from any cause will be encoded as 
modulation of the ZCT data values These data values represent the fluctuations in the 
zero crossing times, and in an ideal system at constant speed, would all be equal to zero. 
Hence, only the fluctuations in load are encoded. Because there are six zero current 
5 crossings in each supply cycle in a typical three phase system, the ZCT sampling is fixed 
at only six samples per cycle. Hence the sampling frequency of the ZCT signal is 300 Hz 
with a 50Hz supply, giving a range of frequencies covering only (0 - 150) Hz in a ZCT 
spectrum. 

The S VAF methodology according to the present invention is not so limited in 
1 0 sampling rate, and hence aliasing may be overcome or minimized, which heretofore has 
been problematic in association with the ZCT technique. Thus, the invention provides 
for deriving sampled data from the rotational motion of the space vector representing the 
three phase currents, instead of just from the zero crossings. The balanced stator currents 
of amplitude I, phase angle (f> and angular frequency a> s for a symmetrical three phase 
1 5 winding are represented by the following equation (2): 
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These stator currents define a space vector as the sum of space vectors of individual 
20 phases given by equation (3): 

h =f |f«(0 + «6(0 + fl 2 'c(0], (3) 

where a = e y2W3 is a space operator. 

25 

For a balanced, symmetrical, steady state system, this vector rotates in the space 
vector plane with constant amplitude in the positive direction at synchronous speed. Its 
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locus in the space vector plane is a perfect circle, as shown in (4): 
i=lei+>eJ"J. (4) 



5 However, for a system, with presence of any other additional component caused 

by either stator fault, rotor fault or other unbalance, of amplitude I C om P .? angular frequency 
6>com P . and phase angle <f> C om P ., the three phase currents may be written as in (5). 
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Substituting (5) into (3), the currents in (5) can be expressed in space vector form: 

(6) 

(7) 



i s = le^e^ + I comp e j ^ e^ 1 



l s Je ^ 1 comp. 



e JC ° compt =1 + 1 



comp. 



The resultant space vector is the sum of the fundamental and additional 
component corresponding space vectors. The locus of the space vector will no longer be 
a perfect circle due to the combined effect of the present components. Apart from 
distortion in the space vector's amplitude, the result is that the rotational speed of the 
space vector fluctuates so that it does not pass any specific point on the circle at equal 
intervals 2^, but does so with fluctuating delays. The space vector angular fluctuation 
(SVAF) method explores these fluctuations by measuring them with respect to the 
balanced referent signal of amplitude / and angular frequency co s . The resultant space 
vector from (6) divided by referent space vector is given in (8): 



(8) 



with /' = Ico ™ p - , and co } =co s -~ co comp . Finally, the angular fluctuation of the space vector 



30 may be found using equation (9): 
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6 = arctan 



f t \ 



f Fsmco't 
= arctan — * 

^ l + /'cosfi)7j 



(9) 



5 The following polynomial expansion (10) computes the arctangent of a variable x, 

when (x < 1): 



arctanW = 0.318253a + 0003314Jc 2 -0.130908x 3 + 0.0068542jt 4 ~0009159;c 5 (10) 
arctan(*)= a • x + p • x 2 - y • jc 3 + 8 • jc 4 - *r • jc 5 



10 Since 7' « 7, this expansion can be applied to (9). If only the first term of the 

expansion, ooc 9 is taken into consideration, together with binomial expansion for the 
argument of the arctangent function in (9), the following expression (1 1) is obtained: 
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(ii) 



Equation (11) indicates frequencies that can be expected in the spectrum of the 
SVAF signal. To obtain the fluctuations in a real system, the instantaneous rotational 
position d x of the current space vector, which advances with each sample instant, is 
measured in the space vector plane. When the expected increment in rotational position 
20 ®[= (o s ti for uniform rotation at angular frequency co s is subtracted, fluctuations 56i in the 
angular position are obtained and given by (12): 

Se^di-mA, (12) 

25 where 6% is the sampled phase angle of the rotating space vector of the real system 

during faulty condition, co s is the supply angular frequency and U is the sample time. A 
voltage signal can be taken as a referent signal, since it will hold accurate information 
about mains supply frequency currently present in the signal. Alternatively, equally 
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spaced samples from the space vector circle can be taken and from two consecutive 
samples and the expected rotation increment 0 can be subtracted, as in equation (13): 

^i=e t -0 t .i-Q. ( 13 ) 

5 

0 is the expected, undisturbed rotation of the space vector between successive 
samples of the space vector taken at intervals t s when the space vector is sampled, where 
0. equals co s t s . Fig. 26 illustrates an exemplary plot 600 of three phase currents 602, 604, 
and 606 with harmonics plotted as current 608 versus time 6 1 0 in the time domain. 
1 0 Instants of time at which samples for SVAF analysis are taken, are indicated with crosses 
612, wherein there are 24 samples taken in the single cycle {e.g. , at 50 Hz) illustrated. 

Fig. 27 illustrates an exemplary plot 620 of Iq component 622 versus Id 
component 624, wherein a circular space vector 626 is illustrated in dashed line for the 
case of an ideal balanced system without harmonics, and where an exemplary space 
1 5 vector 628 is illustrated for the case of three phase system with unbalanced harmonics 
(solid). Samples chosen for the SVAF calculations from the space vector circle are 
marked with crosses 630. Referring also to Fig. 28, an exemplary plot 640 of space 
vector angular fluctuations SVAF 642 versus time 644 illustrates sampled values 5% 646 
presented in the time domain. With the SVAF techniques according to the present 
20 invention, it is possible to select any sampling frequency to sample the space vector 
circle, whereas the sample rate of the zero crossing times ZCT technique is limited. 
Thus, the frequency range in the spectrum of the SVAF signal may be extended, whereby 
aliasing from higher frequencies may be avoided. 

As illustrated in Fig. 29, an exemplary plot 650 of SVAF 652 versus frequency 
25 654 shows an exemplary frequency spectrum 656 of the space vector angular fluctuation. 
The inventors have found that one or more fault indicative frequencies are present in the 
SVAF spectrum 656. In addition, the spectral data includes lines at frequencies 
corresponding to the physical process. For example, torque fluctuations at shaft 
frequency will appear at exactly rotor frequency {f r ), compared to spectra from direct 
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current sampling which provide this information as sidebands to the mains frequency, or 
its harmonics. Other frequencies of interest are illustrated and described in greater detail 
hereinafter. 

Experimental data has been obtained using three induction motors, with special 
5 windings on the stator to enable stator fault experiments. A 4 hp motor was used in 
which slot windings were brought out to panel terminals, to allow short circuiting of 
selected slots to simulate a stator fault. The motor was connected as a four pole induction 
motor with stator windings in 46 slots, each slot having 64 turns. The rotor of the 4 hp 
motor was a wound rotor, which enabled the rotor circuit to be unbalanced by adding 
1 0 extra resistance to simulate rotor faults. Also used were a 2 hp, four pole induction 
motor, and a 2.2 kW, four pole induction motor, in which stator faults were simulated 
through taps on the winding, that could be connected to short circuit two, three, four or 
ten neighboring turns. A resistor was added to the shorting link to limit the fault current 
to protect the winding. In addition, rotor faults were simulated by cutting rotor bars. 
1 5 With respect to stator fault detection, short-circuited turns on the stator of an 

induction motor cause asymmetry of the three phase stator winding. When explained by 
symmetrical components theory, the overall effect of such a fault is the presence of three 
phase negative sequence currents, as is known. Negative sequence currents rotate at 
angular frequency (-o> s ) and so does their corresponding space vector. Thus, the three 
20 phase system from equation (5) has negative sequence currents as additional component 
and angular fluctuations of the resultant space vector from equation (1 1) include terms 
with following frequency 0)'= o s . ( . co^ = 2co s . The spectrum of the SVAF signal will 
include a spectral component at 2f s , where f s is the frequency of power applied to the 
motor. The amplitude of this spectral component has been found by the inventors to 
25 change when a stator fault condition exists. Consequently, the amplitude of the SVAF at 
this frequency may be employed as a diagnostic measure for such stator faults. 

To monitor stator faults in real time, the changes of the 2f s amplitude may be 
monitored. The Goertzel algorithm may be employed to extract this spectral component 
amplitude information using data sampled over a single supply cycle. The Goertzel 
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algorithm advantageously provides expression of the computation of the DFT as a linear 
filtering operation. This may be used, for example, when the DFT is to be computed 
from a sequence of N samples, at a relatively small number M of values, where 
M < log 2 N . In the present example, M= 1 since only one DFT value is needed from the 

5 block of input data. The length of the input data isN=6 for the ZCT method, or more 
for SVAF method, depending on the selected sampling rate. The Goertzel algorithm may 
thus be more efficient than the FFT algorithm for this type of calculation. Where the 
spectral information is thus obtained, e.g., through monitoring the behavior of one 
spectral frequency or component over a period of one mains cycle, a very fast response 

1 0 time results, making it suitable for real time, on-line detection of stator faults. 

An induction motor with 652 turns in each stator phase was simulated. One, two 
three, four and five turns respectively were short-circuited in a computer model and the 
results are illustrated in Fig. 30 as an exemplary plot 660 of 2f s component amplitude 662 
versus time 664. The 2f s diagnostic metric 666 was extracted on-line and its amplitude 

1 5 plotted, wherein the SVAF technique responds virtually immediately to various stator 
faults 670, 672, 674, 676, and 678, corresponding to 1, 2, 3, 4, and 5 of 652 turns shorted, 
respectively, in the motor. The SVAF methodology detects the first supply cycle having 
increased negative sequence current, which is reflected in the diagnostic index 666 (e.g. , 
SVAF spectral component amplitude fluctuation), as illustrated in Fig. 30. The number 

20 of shorted turns, e.g. , the severity of the fault, is also indicated with rise in the amplitude 
of this index 666, wherein the fault was alternatively switched on and off repetitively in 
Figs. 30 and 32. 

Fig. 31 illustrates an exemplary plot 680 of SVAF 682 versus frequency 684, 
centered around the 2f s frequency of 100 HZ, wherein a spectrum 686 is illustrated in 
25 dashed line for a healthy motor, and an SVAF spectrum 688 is illustrated for the motor 
with a stator short circuit (solid), for 1.4% of phase voltage short-circuited turns, as 
obtained experimentally. As illustrated, the amplitude 690 of the component at 
approximately 2fs for the faulted motor is higher than the corresponding amplitude 692 
for the healthy motor. Further experimental results are illustrated in Fig. 32 as an 
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exemplary plot 700 of 2f s component amplitude 702 versus time 704, wherein the 2f s 
diagnostic metric 706 was extracted on-line and its amplitude plotted. The SVAF 
technique responds quickly to various stator faults 708, 710, and 712, corresponding to 
0.7%, 1.05%, and 1.4% of stator turns shorted, respectively, in the motor. 

5 The SVAF technique may further be employed to detect problems associated with 

unbalanced supply voltages, and to distinguish such faults from stator faults in 
accordance with another aspect of the invention. When the supply voltage is unbalanced 
and contains negative sequence voltages, negative sequence currents result on the stator 
of the induction motor. The 2f s current component amplitude is directly affected when 

10 the supply voltage is unbalanced. The invention further provides for distinguishing stator 
faults from supply unbalance conditions based on angular fluctuations in current and 
voltage space vectors. 

For instance, analysis of diagnostic index 2f s from voltage space vector angular 
fluctuations may be combined with analysis of diagnostic index 2f s from the current space 

1 5 vector angular fluctuations, to classify the condition affecting the motor. In this regard, 
the inventors have found that a sudden change in the current 2f s spectral component 
amplitude occurring within one supply cycle with a corresponding sudden change in the 
voltage 2f s spectral component amplitude indicates an unbalanced condition, whereas a 
sudden change in the current 2f s spectral component amplitude occurring within one 

20 supply cycle without such a corresponding sudden change in the voltage 2f s spectral 
component amplitude indicates a stator fault. Thus, fuzzy logic systems or other 
techniques may be employed in the diagnostic component in order to distinguish 
unbalanced power and stator fault conditions in accordance with the present invention. 

The invention further provides for detection and/or diagnosis of rotor problems or 

25 faults within an electric motor using space vector angular fluctuation. The existence of a 
rotor cage fault has been found to cause an electrical asymmetry of the rotor circuit of a 
motor. This asymmetry gives rise to a (l-2s)f s spectral component in the stator current, 
wherein the amplitude of these sidebands reflects the extent of the rotor asymmetry. For 
instance, the interaction of the (l-2s)f s harmonic of the motor current with the 
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fundamental air-gap flux produces speed ripple at 2sf s and gives rise to additional motor 
current harmonics at frequencies given by the following equation (14): 

f rb =(l + 2ks)f S9 * = 1,2,3... (H) 

5 

The three phase system of currents from (5), for rotor fault, includes additional 
components at angular frequency co comp . = (l-2s)co s and the space vector angular 
fluctuation from equation (1 1) holds sinusoidal components at o) ,= Q s - o) comp . = 2sco s and 
its multiples. Spectral analysis of the experimental SVAF signals with rotor asymmetry 

1 0 confirms that this component appears in the spectrum of the SVAF signal and may be 
used as a diagnostic index for rotor faults. In addition, spectral components at sidebands 
(1 - s)2f s and (I + s)2fs, appear around 2f s in the SVAF spectrum and their amplitude 
increases with the severity of the rotor asymmetry. The origin of these sidebands can be 
found when current space vector spectrum is analyzed. For instance, if the motor is 

15 supplied from a non-ideal power source, e.g. , one which has negative sequence current on 
the stator due to the negative sequence voltage, negative sidebands -(l-2s)f s and -(l+2s) 
f s are found in the motor current spectrum. Thus, the SVAF spectrum includes co f = co s - 
cocomp- = <*>s - [-(l+/-2s) 0) s ] = (l+/-s)2co s components. Modulation of the 2f s component 
in the SVAF spectrum by these sidebands causes 2sf s ripple in the 2f s component, giving 

20 an indication of the rotor unbalance, when this frequency component is monitored on- 
line. 

Referring now to Figs. 33-44, an exemplary plot 720 is illustrated in Fig. 33 as 
SVAF power 722 versus frequency 724, wherein a spectrum 726 for a mildly unbalanced 
rotor is shown in dashed line, and a spectrum 728 is shown for a rotor having greater 
25 imbalance, obtained by simulation. For the simulation, different amounts of rotor 

unbalance were created with respect to a healthy rotor resistance of R r '=0.8 1 6 Ohm. The 
first case of unbalance 726 assumes added resistance of AR r '=0.1 Ohm in one rotor 
phase, while in the second case 728, an added resistance of AR r '=0.2 Ohm is used. 
Twice the slip frequency is 6 Hz, and sidebands to the 2f s component occur 
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correspondingly at 94 Hz and 106 Hz for a source frequency of 50 Hz. Fig. 34 illustrates 
an exemplary plot 730 of SVAF power 732 versus frequency 734, wherein a rise in the 
2sf s and 4sf s spectral components is seen due to one broken bar. A normal spectrum 736 
is illustrated in Fig. 34 in dashed line together with a faulted motor spectrum 738. As 
5 illustrated in Fig. 34, the peak 740 for the normal motor is lower than the peak 742 for the 
faulted motor case. 

Fig. 35 illustrates an exemplary plot 750 of SVAF power 752 versus frequency 
754, wherein a rise in the (1-s) 2f s and (1+s) 2f s sidebands is illustrated from the mild 
rotor resistance imbalance spectrum 756 to a spectrum 758 for higher degree of 

10 imbalance, for the imbalance conditions discussed above with respect to Fig. 33. In Fig. 
36, a plot 760 of SVAF power 762 versus time 764 illustrates oscillations of the 2f s 
component at 2sf s frequency, via mild imbalance curve 766 and higher imbalance curve 
768 for the same degrees of rotor resistance unbalance. Fig. 37 provides a plot 770 of 
current space vector power 772 versus frequency 774 illustrating an exemplary current 

15 space vector spectrum 776, and Fig. 38 provides a plot 780 of voltage space vector power 
782 versus frequency 784 illustrating a voltage space vector spectrum 786. 

Figs. 37 and 38 illustrate the spectrum of the current and voltage space vectors, 
respectively, under unbalanced rotor (AR/=0.1 Ohm) and unbalanced voltage supply 
condition of 0.6%, wherein supply voltage unbalance condition is the ratio of negative to 

20 positive sequence voltage. This voltage negative sequence causes 4.86% unbalance in 
the stator currents. A spectral component at (-50 Hz) is present in the spectrum of the 
voltage space vector, as well as in the spectrum of the current space vector. Negative 
sidebands -(l-2s) f s and -(l+2s) f s in the current spectrum are also distinguishable. 
These cause sidebands (1-s) 2f s and (1+s) 2f s in the SVAF spectrum, which are shown on 

25 the Fig. 39, which provides a plot 800 of SVAF power 802 versus frequency 804. As 
illustrated in the spectrum 806 of Fig. 39, spectral components 8 10, 812, 814, 816, 818, 
820, 822, 824, 826, 828, and 830 are seen for 2sf S9 4sf s , 6sf S9 (l-s)6£, (l-s)4/„ (l-$)2f S9 2f S9 
(l+s)2f s , (l+s)4/i, and (l+s)6/i, respectively. 

Fig. 40 illustrates a plot 840 of SVAF power 842 versus frequency 844 for a 
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normal motor spectrum 846 and a broken bar fault spectrum 848, wherein fault indicative 
frequency spectrum components 850 are illustrated. Results presented are for normal 
condition on the motor and for one broken rotor bar 5 out of twenty six. Supply voltage is 
1% unbalanced. Figs. 34 and 40 show 2sf S9 and (l-s) 2f s and (1+s) 2f s spectral 
5 components, respectively. Fig. 41 illustrates a plot 860 of 2f s component amplitude 862 
versus time 864 for a normal motor 866 and for a motor having one broken bar 868. 
Finally, complex spectra of the current and voltage space vectors are presented on Figs. 
42 and 43, respectively. Fig. 42 includes a plot 870 of current space vector power 872 
versus frequency 874 showing a spectrum 876, and Fig. 43 illustrates a plot 880 of 

10 voltage space vector 882 versus frequency 884 with a spectrum 886. Negative sequence 
current and voltage, as well as -(l-2s)f s negative current sideband can be clearly be seen. 
Fig. 44 presents full spectrum with rotor fault indicative frequencies, wherein a plot 890 
of SVAF power 892 versus frequency 894 illustrates a spectrum 896. 

The angular fluctuations of the induction motors current space vector may thus be 

15 employed in accordance with the invention as diagnostic data for stator and rotor 
induction motor faults, as well as for diagnosing unbalanced supply conditions. The 
SVAF methodology allows monitoring of fault indicative changes in the diagnostic data, 
through analyzing characteristic frequencies in the spectrum of the space vector angular 
fluctuations. The following table summarizes various diagnostic indices for stator and 

20 rotor faults which may be detected according to the invention: 



TABLE I 

DIAGNOSTIC INDICES FOR SVAF METHOD 



Diagnostic 
index 


Origin of 
index 


Require 
d 

for 


Monitoring type 


Trend 


Current 


Current 


Stator 


Real time 


Increase with neg. 


SVAF 2f s 


neg. 


fault 




sequence 




sequence 








Voltage 
SVAF 2f s 


Voltage 
neg. 

sequence 


Stator 
fault 


Real time 


Increase with neg. 
sequence 
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Current 


(l-2s)f s 


Rotor 


On-line spectrum 


Increase with 


SVAF 2s£ 


current 


fault 




fault 




siaeoana 








Current 


-(l-2s)f s 


Rotor 


On-line spectrum 


Increase with 


SVAF 


current 


fault 




fault 


(l-s)2f s 


sideband 








Current 


-0+2s)f s 


Rotor 


On-line spectrum 


Increase with 


SVAF 


current 


fault 




fault 


(l+s)2f s 


sideband 









Stator faults may thus be indicated in the spectrum as change of the 2f s and other 
component amplitudes. The Goertzel algorithm may be advantageously employed to 
extract the component's amplitude from the signal during each cycle of the mains supply. 
5 This avoids lengthy procedure required for creating frequency spectrum by Fourier 
Transformation as well as allowing real-time monitoring of motor faults. The 2f s 
component will respond to the change in the angular fluctuations as soon as fault 
happens, allowing fast recognition in accordance with the present invention. The 
recognition process may include angular fluctuations of the voltage space vector, since 

10 voltage unbalance and load changes will affect this frequency. Rotor faults are indicated 
in the current space vector angular fluctuations spectrum with spectral lines at 2sf s caused 
by (1 - 2s)f s and (1 + 2s)f s in the motor current spectrum, and at (1 - s)2f s and (1 + s)2f S9 
which are directly caused by spectral lines ~(1 - 2s)f s and -(1 +2s)f s in the motor current 
space vector spectrum. Additionally, the 2f s SVAF spectral component will oscillate at 

15 2sf s when monitored during each cycle, giving indication of the rotor unbalance in real- 
time. 

Although the invention has been shown and described with respect to certain 
illustrated aspects, it will be appreciated that equivalent alterations and modifications will 
occur to others skilled in the art upon the reading and understanding of this specification 
20 and the annexed drawings. In particular regard to the various functions performed by the 
above described components (assemblies, devices, circuits, systems, etc.), the terms 
(including a reference to a "means' 1 ) used to describe such components are intended to 
correspond, unless otherwise indicated, to any component which performs the specified 
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function of the described component (e.g., that is functionally equivalent), even though 
not structurally equivalent to the disclosed structure, which performs the function in the 
herein illustrated exemplary aspects of the invention. In this regard, it will also be 
recognized that the invention includes a system as well as a computer-readable medium 

5 having computer-executable instructions for performing the acts and/or events of the 
various methods of the invention. 

In addition, while a particular feature of the invention may have been disclosed 
with respect to only one of several implementations, such feature may be combined with 
one or more other features of the other implementations as may be desired and 

10 advantageous for any given or particular application. As used in this application, the 
term "component" is intended to refer to a computer-related entity, either hardware, a 
combination of hardware and software, software, or software in execution. For example, 
a component may be, but is not limited to, a process running on a processor, a processor, 
an object, an executable, a thread of execution, a program, and a computer. Furthermore, 

15 to the extent that the terms "includes", "including", "has", "having", and variants thereof 
are used in either the detailed description or the claims, these terms are intended to be 
inclusive in a manner similar to the term "comprising." 
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